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PERIODIC PERTURBATIONS OF CONSTRAINED MOTION 
PROBLEMS ON A CLASS OF IMPLICITLY DEFINED 

MANIFOLDS 

ALESSANDRO CALAMAI AND MARCO SPADINI 



Abstract. We study forced oscillations on differentiablc manifolds which are 
globally defined as the zero set of appropriate smooth maps in some Euclidean 
spaces. Given a T-periodic perturbative forcing field, we consider the two 
different scenarios of a nontrivial unperturbed force field and of perturbation of 
»"H ' the zero field. We provide simple, degree-theoretic conditions for the existence 

of branches of T-periodic solutions. We apply our construction to a class of 
■^T , second order Differential-Algebraic Equations. 
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1. Introduction 



In this paper we study T-periodic solutions of some parametrized families of 
T-periodic constrained second order Ordinary Differential Equations (ODEs). In 

rS^ | other words, we study forced oscillations on a diffcrcntiable submanifold of some 

Euclidean space. As a physical interpretation, the equations we consider here rep- 
resent the motion equations of a constrained system, the manifold being the con- 
straint. We work under the assumption that such a manifold is globally defined 
as the zero set of a C^-smooth map and, given a T-periodic perturbative forcing 
field, we consider two different scenarios according whether there exists a nontrivial 

_j_ \ unperturbed force field. Namely, let g: U C R m x R s — > W be a C°° map such 

that the Jacobian matrix of g with respect to the last s variables, d2g(x, y), is non- 
singular for all (x,y) E U. This implies that € R s is a regular value of g so that 
M = g _1 (0) is a C°°-smooth submanifold of R k = W n x W of dimension m (and 

■^- \ codimension s). 

For second order ODEs on differentiable manifolds we adopt the notation of e.g. 
0]. We consider parametrized second order ODEs on M that, with this notation, 
assume the following forms: 



(1.1a) 6r = /K,0 + Afc(t,£,0 

and 



X 

.?.' (Lib) ir = M(t,t,£), 

where A > is a parameter, £(£) £ M, £„■ stands for the tangential part of the 
acceleration, h: R x TM —} R fc and /: TM — > R fc are continuous maps with the 
property that /(^,r?) and h(t,£,r)) belong to T^M for any (t,£,rj) E R x TM, and 
h is T-periodic in the first variable. 

Notice that the assumptions on M and g imply that, locally, the manifold M 
can be represented as graph of some map from an open subset of R m to R s . Thus, 
writing £(t) = (x(t),y(i)), Equations (|1.1[) can be locally simplified via the implicit 
function theorem. In view of this fact one might think that it is possible to reduce 
Equations (|1.1|) to ordinary differential equations in R m . It is not so. In fact, 
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2 ALESSANDRO CALAMAI AND MARCO SPADINI 

globally, M may not be the graph of a map from an open subset of R m to R s as, 
for instance, when U = R 3 and g : R x R 2 — > R 2 is given by 

9{x, y) = g{x; j/i, y-x) = (e yi cos(y 2 ) - a:, e yi sin(y 2 ) - a;) . 

In this case, although det d2g(x,y) ^ 0, one clearly has that the 1-dimensional 
manifold M = g~ l {Qi) is not the graph of a function x n- (yi(x), yi{xy\ . In fact, M 
consists of infinitely many connected components each lying in a plane y 2 = j + lir 
for I e Z. 

Observe also that even when M is a (global) graph of some map T, the expression 
of r might be too complicated to use or impossible to determine analytically so 
that, the simplified versions of Equations (ll.lj) may be too difficult to use. A 
simple example of this fact is obtained by taking m = s = l,[/ = RxR and 
g(x,y) = y 7 + y-x 2 + x 5 . 

A pair (A, £), with A > and £: R —> M a T-periodic solution of (jl.lal) (resp., of 
(|l.lb|l ) corresponding to A, is called a solution pair of (jTTTaJ) (resp., of (jl.lb|l ). The 
set of solution pairs is regarded as a subset of [0, +oo) x C^(M), where C^(M) is 
the metric subspace of the Banach space C^(R fe ) of the T-periodic C 1 maps from 
R to M. In this paper we investigate the structure of the set of solution pairs of 
Equations (|l.lj) . Namely, we will prove global continuation results for solution pairs 

of jnj. 

Our results here stem from a combination of those of [SJ [TU] about branches 
of T-periodic solutions of constrained second order equations, with the formula 
contained in [5] for the computation of the degree of a tangent vector field in terms 
of the Brouwer degree of an appropriate map. In the situation considered in this 
paper (namely, for manifolds which are of the form M = flf (0) with g as above), 
our results are complementary and, in some sense, improve upon those of [SJ 110) 
since the Brouwer degree of the considered maps, having a simpler nature than the 
degree of a tangent vector field, is in principle easier to compute. In fact, the degree 
of a tangent vector field can be seen as an extension of the notion of Brouwer degree 
(a different one being the notion of degree of maps between oriented manifolds). 
Actually, since vector fields in Euclidean spaces can be regarded as maps and vice 
versa, the degree of a vector field is essentially the Brouwer degree, with respect to 
0, of the field seen as a map. 

A remarkable application of our results, which justifies our interest in the par- 
ticular differentiable manifolds considered in this paper, is the content of the last 
section where we study periodic perturbations of a particular class of second order 
Differential-Algebraic Equations (DAEs) in semi-explicit form. 

Recently, DAEs have received increasing interest due, in particular, to applica- 
tions in engineering and have been the subject of extensive study (see e.g. [T3] for 
a comprehensive treatment) aimed mostly (but not only) to numerical methods. 
Our approach here, as in our recent works [H [^J [TS], is directed towards qualitative 
theory of some particular DAEs which are studied by means of topological meth- 
ods, making use of the equivalence of the given equations and suitable ODEs on 
manifolds. In particular, in this paper we investigate particular DAEs of second 
order, whereas the papers [TJ ® HH] were devoted to the first order case. The type 
of equations that we study here may be used to represent some nontrivial physical 
systems, as, for instance, constrained systems (see e.g. .17 ). As we will prove, such 
equations are equivalent to second order ODEs on precisely the kind a manifold that 
we consider in the first part of the paper. This fact enables us to take advantage 
of the continuation results in [SJ [TU] and to find some surprisingly simple formulas 
useful for the study of connected sets of forced oscillations when a periodic forcing 
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term is introduced. Some applications of our results are shown in a few illustrative 
examples. 

We observe that our results, besides their intrinsic interest, can be useful in 
order to establish existence theorems in presence of a priori bounds (for a general 
reference on continuation methods, see e.g. |16|). In addition, they can be used to 
get some 'topological' multiplicity theorems for forced oscillations as in [7]. We will 
not pursue these lines here, though. 

Finally, we wish to point out that in this paper we make use of deep results which 
are based on the fixed point index, but our techniques require just the notion of 
the well-known Brouwer degree. Our effort has been to make the paper accessible 
also to researchers which are not particularly familiar with topological methods. 

2. Preliminaries 

2.1. Tangent vector fields and the notion of degree. We now recall some 
basic notions about tangent vector fields on manifolds as well as the notion of 
degree of an admissible tangent vector field. 

Let M C M. k be a manifold. Let w be a tangent vector field on M, that is, 
a continuous map w: M — > M. k with the property that w(£) £ T^M for any £ G 
M. If w is (Frechet) differentiable at £ £ M and w(£) = 0, then the differential 
dw^: T^M — > R fe maps T^M into itself (see e.g. [15]), so that the determinant 
det du>£ of du>£ is defined. If, in addition, £ is a nondegenerate zero (i.e. dw^ : T^M — > 
R k is injective) then £ is an isolated zero and det dw^ ^ 0. 

Let W be an open subset of M in which we assume w admissible (for the degree); 
that is, the set u> _1 (0) P\W is compact. Then, one can associate to the pair (u>, W) 
an integer, deg(w,W), called the degree (or characteristic) of the vector field w in 
W, which, in a sense, counts (algebraically) the zeros of w in W (see e.g. [5l fT5lfT5] 
and references therein). In fact, when the zeros of w are all nondegenerate, then 
the set iu -1 (0) n W is finite and 

(2.1) deg(ui, W) = >J signdetdtuj. 

fe«;- 1 (o)nw 

Observe that in the flat case, i.e. when M = M. k , deg(w, W) is just the classical 
Brouwer degree with respect to zero, deg B (w,V,0), where V is any bounded open 
neighborhood of w _1 (0) n W whose closure is contained in W. 

The notion of degree of an admissible tangent vector field plays a crucial role 
throughout this paper. All the standard properties of the Brouwer degree for con- 
tinuous maps on open subsets of Euclidean spaces, such as homotopy invariance, 
excision, additivity, existence, still hold in this more general context (see e.g. [5]). 

Remark 2.1. The Excision Property allows the introduction of the notion of index 
of an isolated zero of a tangent vector field. Indeed, let £ £ M be an isolated 
zero of w. Clearly, deg(u>, V) is well defined for each open V C M such that 
V n iu _1 (0) = {£}. By the Excision Property deg(w, V) is constant with respect 
to such V 's. This common value of deg(w, V) is, by definition, the index of w at 
£, and is denoted by i(u>,£). Using this notation, if (w, W) is admissible, by the 
Additivity Property we get that, if all the zeros in W of w are isolated, then 

(2.2) deg(w,W)= J2 i ( w >£)- 

5eu)- 1 (o)nw 

By formula (|2.1|) we have that, if £ is a nondegenerate zero of w, then 

i (w, £) = sign det dw^ . 
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Notice that (|2.1[) and (|2.2p differ in the fact that, in the latter, the zeros of w are 
not necessarily nondegenerate as they have to be in the former. In fact, in (J2.2D , w 
need not be differ entiable at its zeros. 

2.2. Tangent vector fields on implicitly defined manifolds. Let \1/: MxM-> 

R fe be a (time-dependent) tangent vector field on M C R fe that is, a continuous 
map with the property that ^(i, £) G T^M for each (£,£) € R x Af. Assume that 
there is a connected open subset U of R fc and a C°° map <? : U — ?> R s with the 
property that M = g (0). Suppose that up to change of coordinates, writing 
R fc = R m x R s , m = k — s, one has that the partial derivative of g with respect to 
the second variable, d2g(x,y), is invertible for each (x,y) G f . 

According to the above decomposition of M fc , we can write, for any £ G R fc , 
£ = (je, y) and, with a slight abuse of notation, for any t G R 

*(t,0 = *(t,x,y) = (*i(t,s,i/),* 2 (*,a;,»)). 

Notice that one must have 

(2.3) *a(*,ac,i/) = -(ftsfol/jr^sfoy^ifcaM/). 

In fact, ^J(£, £) G T^M being equivalent to ^(t, £) G ker g'(x,y), one has for each 
(t, x, y) G R x M that 

= #'(*, j/)*(t, *, y) = 9isr(«, y)*i(i, *> y) + d 2 g(x, y)* 2 (t, *, y), 

which implies (|2.3|l : here g'(x,y) denotes the Frechet differential of g at (x,y). 

We now give a formula for the degree of tangents vector fields on M in terms of 
(potentially easier to compute) degree of appropriate vector fields on U. 

Remark 2.2. Assume that tp: M — > R fc is a tangent vector field on M. Since 
M = g _1 (0) is a closed subset of the metric space U, the well-known Tietze's 
Theorem (see e.g. [3],) implies that there exists an extension ip: U —¥ M fe ofip. 



Remark 12.21 shows that it is not restrictive to assume, as we sometimes do, that 
the given tangent vector fields are actually defined on a convenient neighborhood 
of the manifold M. In fact, although an arbitrary extension of "0 may have many 
zeros outside M, we are interested in the degree of tp on M which only takes into 
account those zeros of ip that lie on M. 

The following result (see [2], [18]) gives a formula for the degree of tangents 
vector fields on M in terms of degree of appropriate vector fields on U. 

Theorem 2.3. Let U C R m x R s be open and connected, let g: U -^ W be a C°° 
function such that d2g(x,y) is nonsingular for any (x,y) G U and let M — g -1 (0). 
Assume that tp: M — > R m x R s is tangent to M and let tpi be the projection on 
R fc of an arbitrary continuous extension tp of tp to U. Define T: U — » R m x R s by 
F{x,y) = (tpi{x, y), g(x, y)) . Then, JF is admissible in U if and only if so is tp in 
M, and 

(2.4) deg(^,Af)=sdeg(^,(7), 

where s is the constant sign o/det d2g(x,y) for all {x,y) G U . 

3. Second order equations on manifolds 

As in the previous section, let M be a submanifold of M fc . Given a continuous 
map ip: R x Till ->■ M fe such that tp(t,^Tj) G T^M for any (t,£,r]) G R x TM, we 
will say that tp is tangent to M. 

It is known (compare [4]) that the motion equation associated with the force (p 
can be written in the form 

(3.1) S« = <p(t,£,£), 
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where £„• stands for the parallel component of the acceleration £ £ R fc at the point 
£. Namely, £„. denotes the orthogonal projection of £ onto TjAf . A solution of 
(|3.ip is a C 2 map £: J — > M, defined on a nontrivial interval J, such that £*■(*) = 
¥>{t,£,(t),£(t)) for alHe J. 

Further, Equation (|3.1[) can be written in an equivalent way as a first order 
equation on the tangent bundle TM in the form 

(3.2) C = *(*,0. 
where, for £ = (£, 77), 

(3.3) *(*;£,»?) = (»?, K£,»7)+¥>(t; &»?)), 

and the map r: TAf — ► R fc is smooth, quadratic in the second variable f] £ TjAf 
for any £ £ Af, and with values in (T^M)- 1 . Such a map r is strictly related to 
the second fundamental form on M and may be interpreted as the reactive force 
due to the constraint M. Actually r(£, 77) is the unique vector in R fc which makes 
{jli r (£i r ?)) tangent to TAf at (£,77). It is well known that \1/, called the second 
order vector field associated to (p, is a (time-dependent) tangent vector field on 
TM. Hence (|3. 21) is actually a first order equation on TM. 

We now focus on the computation of the map r. 

Let £ : J C R -> M be a local solution of (13. ip taking values in a neighborhood 
V C M of £0 := £(£0) f° r some to £ J. Restricting V if necessary, we can assume 
that there exists U C R fc and a C°° function G: 17 -* M s such that (7 n M = V = 
G _1 (0). Differentiating twice at £0 the relation G(£(£)) = 0, we get 

(3.4) G'%)(£ ,eo)+G'(£o)|(£ )=0, 

where, for the sake of simplicity, we write £0 instead of £(to). According to Equa- 
tions (I^2l)-(l3~3l) 



(3.5) £(to) = ¥>(to,6>,&)+r(eo,&)- 

Since £,(*„) = </?(£o,£o, £0) G I&M, one has G'(£ )¥>(£o,£o,£o) = 0. Hence (0 
yields 

G'(eo)r(eo,eo)=-G"(Co)(eo,eo). 
And, finally, 

(3.6) r(£ ,£o) = - (G'(eo)|(T, A/)-) _1 G"(eo)(eo,eo)- 

Recall, in fact, that G'(£o)|(T e Af)- 1 - is invertible because Tf Af is the kernel of G'(£o)- 
Observe that equations (|3.1[> (|3.5p have another immediate interesting conse- 
quence. Namely, that if 1 1— > £(£) is an Af -valued G 2 curve then 

(3-7) ? > c(t)C(t)=r(C(t),C(*)) ) 

where / P^7 t - ) denotes the orthogonal projection of R fc onto the orthogonal comple- 
ment to Tf( t )Af. 

When Af is of the form considered in Sect ion |2~21 we are able to get a more explicit 
expression for the inertial reaction r. The expression that we are going to obtain, 
although strictly speaking not necessary for the carrying out of our arguments, 
can be useful for understanding our geometric setting. The reader may skip the 
remaining part of this section at a first reading. 

Assume, as in Section [2~2l that M = <7 _1 (0) where g: U — > R s is a G°° map with 
d2g(x, y) is invertible for each (x, y) £ U. Let £ : J C R — > Af be a local solution of 
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(|3.1[) and write £ = {x,y). Setting £o = (xo>2/o) and Co = (xo,yo)> Equation (|3.4[> 
becomes 

(3.8) ff"(xo,yo)((*o,S'o))(*o,yo)) + dig(x ,yo)x(t Q ) + d 2 g{xo,yo)y(t ) = 0. 

Thus, r(xoj2/o; ioj^o) is the (unique) solution lying in (T( Xo y )M) of the linear 
equation 

(3.9) g'{x Q ,y Q )r{x a ,y Q ;x ,y ) = -g"(x , yo)((x , yo), (xq, yo))- 

Now, to find the explicit expression of r we need two elementary linear algebra 
lemmas. The proof of the first is straightforward and is left to the reader. 

Lemma 3.1. Let A be an s x m matrix and B be an s x s matrix. Let W — 
{(u,v) € R m xR s : Au + Bv = 0}. Assume B invertible. Then, W 1 - = {(u,v) <E 
l m xl s :u = A T (B- 1 ) T v}. 

Lemma 3.2. Let A be an s x m matrix and B be an s x s matrix. Assume B 
invertible. Then, the s x s matrix C = AA T (B~ l ) T + B is invertible too. 

Proof. To prove that C is invertible we will show that so is B~ X C . To see this, we 
prove that B~ l C is symmetric and positive definite. Notice that 

B~ X C = B- 1 AA T (B- 1 ) T + 1 = B- 1 A{B- 1 A) T + /, 

thus B~ X C is of the form XX T + L with X a s x m matrix, / being the identity. 
In particular, B~ X C is symmetric. Now, observe that XX T is always symmetric 
and positive semidefinite and thus B~ X C is positive definite, being the sum of a 
positive definite matrix (the identity) and a positive semidefinite matrix. □ 

Set now A — dig(xo,yo), B = d2g{xo,yo), and C = AA T (B^ 1 ) T + B (recall 
that B is an invertible s x s matrix by assumption). Further, let r = (u,v) and 
a = —g"(xo,yo)((%o>yo)> {xo,yo)), so that Equation (|3.9[) can be written as 

(3.10) Au + Bv = a. 

Taking into account the above lemmas, simple calculations show that Equation 
(|3.10p has a unique solution (u, v) 6 W n x R s that lies in (Tr Xg Vo )Mj and is given 

by 

u = A T (B- 1 ) T C- 1 a, 

v = C- 1 a. 

Thus, r(xo,yo]&o,yo) has the following expression: 

(3.11) r{xQ,y Q ;x ,y ) = I c _ x ^ 

where a = -g"(x ,yo)((xo,yo), (xo,yo))- 

Example 3.3. Let k = 2, m = 1 and g(x,y) = tj^ 2 — y — 2. We wish to compute 
the map r in this case. Clearly M = <7 _1 (0) is a parabola. With the above notation, 
we have 

A = d ig (x,y)=x, B = d 2 g(x,y) = -1, C = AA T (B- 1 ) T + B = 1 - x 2 . 
Also, 

a = -g"(x iy )((l),(l))=-(*y)^ °)(|)=-± 2 . 
Therefore, 



r(x,y;x,y) 



A T (B~ 1 ) T C- 1 a 



C- x cj 
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Example 3.4. Let k = 3, m = 2 and g{x,y) = z — x 2 — y 2 . We wish to compute the 
map r in this case. Clearly M = 5 _1 (0) is a paraboloid. With the above notation, 
we have 

A = d ig (x, y) = ( -2* -2») , B = d 2 g(x, y) = l, 

\TlT3-l\T 



Also, 



C = AA 1 (B- 1 ) 1 +B = Ax z + Ay z + 1. 



■g"(x,y,z) ((?)'(?)) = -(*»*) 



-2 











-2 















(S) 



= 2x 2 



US 



Therefore, 



r(x,y,z;x,y,z) 



A T (B- 1 ) T C- 1 a 

c-v 



/ -4x(x 2 + y 2 )\ 
1 + 4a; 2 + Ay 2 
-Ay(x 2 + y 2 ) 



1 + Ax 2 + Ay 2 
2x 2 + 2y 2 



The examples above, in conjunction with Equations (13.21) and (|3.3p . show that 
even for simple manifolds second order equations can be rather complicated. The 
following example confirms this fact, strengthening the case for the indirect methods 
investigated in this paper. 

Example 3.5. Let k = 3, m = 1 and g(x, y, z) = (z 3 + z — x. z — y + a; 2 ) . We wish 
to compute the map r in this case. Here M = g~ l (0) is a smooth curve. With the 
above notation, we have 



A = dig{x,y) 



-1 

2x 



B = d 2 g{x,y) 



Also, 



° = -g"(x,y,z) ((?),( 
Lengthy but straightforward computations yield 




-1 



6i 2 
2x 2 



3z 2 + l 
1 



r{x,y,y;x,y,z) 



2(1 + 4a; 2 



20a; 2 + 16a; 4 + 18z 4 + 36z 4 a; 2 + 12z 2 + 24z 2 a; 2 



'(18a;z 4 



- 3z 2 + 12a;z 2 + 8a; 3 + 4a; + l)x 2 - (18a;z 2 + 6a; + 12a; 2 
+9i 2 z 2 + 3i 2 - 2a; 2 - 4i 2 a; 2 - 9a; 2 z 4 - 6£ 2 z 2 
18z 2 z 2 + 6i 2 + 36i 2 z 2 a; 2 + 12i 2 a; 2 + x 2 + 4x 2 x 2 



-9)i 2 



3.1. Branches of harmonic solutions. In what follows, if N is a differentiable 
manifold embedded in some R fe , we will denote by C^(N), the metric subspace of 
the Banach space (C^(R k ) , | • |i) of all the T-periodic C 1 maps £: K ->• N with 
the usual C 1 norm. 

In this section we recall two known results from [5] and [10. about the sets of 
solution pairs of (II. lap and of (|l.lbp . i.e. of those pairs (A, £) € [0, oo) x C^(N) 
with £ a T-periodic solution of (|l.la|) and of (ll.lbp . respectively. Recall that a 
solution pair (A, £) of (ll.lap is said to be trivial if A = and £ is constant, whereas 
a solution pair (A, £) of (jl.lbp is trivial merely if A = 0. 
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For the sake of simplicity we make some conventions. We will regard every space 
as its image in the following diagram of natural inclusions 

N ► C^(N) 



(3.12) 

[0,oo) x N > [0,oo) x C^(N) 

In particular, we will identify N with its image in C^,(N) under the embedding 
which associates to any £0 G N the map £0 G @t(N) constantly equal to £o- 
Moreover we will regard N as the slice {0} x N C [0,oo) x N and, analogously, 
C)p{N) as {0} x Cj.(N). We point out that the images of the above inclusions are 
closed. According to these identifications, if Vt is an open subset of [0, 00) x Cj,(N), 
by N n f2 we mean the open subset of N given by all £0 G N such that the pair 
(0,^o) belongs to fi. If O is an open subset of [0, 00) x N, then OnJV represents 
the open set {£0 G N : (0, Co) G O). 

We need to introduce some further notation. Given / : TN — > R fc tangent to the 
manifold N, we define the tangent vector field /|jv: N — > Mr given by /|jv(£) — 
/(£, 0) for any £ £ N. The following results of [5] and of [10] play a central role: 

Theorem 3.6 ([ID] Th. 4.2). Let N C M*-' fee a boundaryless manifold, and let 
f: TN ->• M fe and /i: K x TiV ->• M* 6e tangent to N. Assume that h is T '-periodic 
in the first variable. Let ft be an open subset of [0, 00) x C^(N). Assume that 
deg (/|jVi N n Cl) is well defined and nonzero. Then, fl contains a connected set T 

of nontrivial solution pairs for (jl.la[) whose closure meets (/|jv) (0) H fi and is 
not contained in any compact subset ofVt. Ln particular, if N is closed in M. k and 
il = [0,oo) x C)p{N), then V is unbounded. 

Given h : M x TN — > R k tangent to the manifold N, we define the mean value 
tangent vector field u>h '■ N — J> M. k given by 

1 - T 



(3.13) Wh ($ = -J h(t,$,Q)dt. 

Theorem 3.7 ([S] Th. 2.2). Let N C M. be a boundaryless manifold, and let 
h: Rx TN — > M. k be tangent to N and T -periodic in the first variable. Let Vt be an 
open subset of [0, 00) x C^(N). Assume that deg (wh,N n f2) is well defined and 
nonzero. Then, $7 contains a connected set T of nontrivial solution pairs for (jl.lbj) 
whose closure meets w^ (0) n £1 and is not contained in any compact subset of Vt. 
Ln particular, if N is closed in M. k and il — [0, 00) x Cj,(N), then T is unbounded. 

4. Main results 

Throughout this section U will be an open and connected subset of R m x K s . 
From now on, we will denote by Pri the projection onto the first factor in the 
Cartesian product l ra xM s . 

We will always assume that g : U —t M. s is a C°° function such that d2g(x,y) is 
nonsingular for any [x, y) G U, and that M = g~ 1 (0). Moreover we will identify 
£ = (x, y) G U. It will also be convenient, given a continuous tangent vector field, 
/ : TM — > M. m xR s , to denote by / an arbitrary extension of the map /|m : (x, y) <—> 
f(x, y, 0, 0) to U (see Remark |2~2^) and to let Pri f(x, y) be the projection of f(x, y) 
on M. m for any (x, y) G U. 

Theorem 4.1. Let f : TM -> R m x R s and h:Rx TM -* M m x W be continuous 
tangent vector fields, with h of a given period T > in the first variable. Define 
T: U -> M m x R s by T(x,y) = (Pi'i f(x,y),g(x,y)) for any £ = (x,y) e U. Given 
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an open set f2 C [0, oo) x C^(M), let O C R m x W be open with the property that 
O PI M = M n fi. Assume that deg(_F, 0) is we// defined and nonzero. Then there 
exists a connected set T of nontrivial solution pairs for (jl.lap in f2 whose closure 
in $7 meets (/Uf) _1 (0,0) l~l and is not compact. In particular, if M is closed in 

R m xR s and Q = [0, oo) x C^{M), then V is unbounded. 

Proof. By Theorem 12.31 we have 

|deg(/|M,MnO)| = |deg(/|M,0nM)| = |deg(^,O)|. 
Thus, dcg(/, M Plfi)^0 and the assertion follows from Theorem 13.61 □ 



The following result concerns Equation (jl.lbl) . As in the previous section, given 
h: M. x TM —> R k tangent to M, we define the tangent vector field Wh on M by 
(|3.13[) . Moreover by Pri Wh(x,y) we denote the projection on R m of an arbitrary 
extension Wh{x,y) of the map Wh{x,y) to U. 

Theorem 4.2. Let ft,: E x TM —> M." 1 xW be a continuous tangent vector field, 
of a given period T > in the first variable. Define $: U — > M. m x R s by &(x, y) — 
(Pri Wh(x, y), g(x, y)) for any £ = (x,y) £ U. Given an open set fi C [0, oo) x 
C\,{M), let O cW n x W be an open subset with the property that OnM = Mnfi. 
Assume that deg($, O) is well defined and nonzero. Then there exists a connected 
set r of nontrivial solution pairs for (jl.lbp in fl whose closure in fl is not compact 
and meets the set w^ (0) n f2. In particular, if Ad is closed in M. m X M. s and 
£1 = [0,oo) x C)p{M), then T is unbounded. 



Proof. By Theorem 12.31 we have 

IdegK.M nn)| = \deg(w h ,OnM)\ = |deg(*,0)|. 

Thus, deg(wh, M n il) ^ and the assertion follows from Theorem 13.71 □ 

Remark 4.3. // is worth noticing that, in Theorem \4- 1\ a point (x,y) £ M. m x M s 
belongs to (/|m) (0,0) if and only if F(x,y) = 0. Similarly, in Theorem \4-ty 
(x, y) £ w^ 1 (0) if and only if $(x. y) = 0. 

Example 4.4. Consider the simple case when m = s=l, (7=1x1 and 
g(x, y) = 3 (x 2 + \)y + ^y 3 + x. Assume that a point P of mass m is constrained 
without friction to the curve M — (7 _1 (0) lying on a vertical plane (see figure\7}) 
and acted upon by the gravitational force ip(x,y) = (0, — mg) T . Here, g denotes 
the gravitational constant. The motion of P is determined by the following second 
order equation on M: 

where (x, y) and V^ denotes the orthogonal projection ofM 2 onto T^M . In order to 
apply Theorems \4-l\ or \4-.2\ it is necessary to compute V^ip(^). This is readily done, 
in fact, 

/(e):= ^ (0 = ^ ) _ MO,vff 

l|Vg(Olr 

where '{•,•}' denotes the scalar product. In coordinates, we get the rather unappeal- 
ing formula below: 

„ s Zrng ((y 2 + 3x 2 + 3)(2xy + 3) 

f(x,y)- ' 



(6xy + 9) 2 + (y 2 + 3x 2 + 3) 2 V -3(2a;y + 3) 
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/(£(*)) 




V / 2 + 2\/2 



rV / 2 + 2 % /2 



gravity 



r %/2 + 2\/2 



\/2+2v^ 



FIGURE 1 . The setting of Example 14.41 



Simple considerations show that there are exactly two zeros of f(x, y) in U cor- 
responding to the two real roots of the polynomial 4x 4 — 4a; 2 — 1, that is x — 
ii-v/2 + 2-v^. Now, letO = (-00,0) xl,!l= [0, oo) x C^(0), and 



F(x,y) 



(3mg{y 2 + 3x 2 + 3){2xy + 3) 1 



(x 2 + l)y + 



1 



y 3 + x 



\(6xy + 9) 2 + (y 2 + 3x 2 + 3) 2 ' 3 v " ~ y " ' 27 £ 
^4 lengthy but straightforward computation shows that deg(J r , O) = 1. Hence, The- 
orem \4-i\ yields that if any T -periodic force Xh(t, £, £), A > 0, is added to f there 
exists an unbounded connected set T of nontrivial solution pairs for (|l.la|) in Q 
whose closure in fi meets / _1 (0) n fl. 




Figure 2. The setting of Example 14.51 



Example 4.5. Consider a unit mass point P constrained without friction to the 
parabola M of Example \3.!A of equation \x 2 — y — 2 = lying on an horizontal 
plane (so that the weight has no effect on the motion). Assume that the extremities 
of a spring with elastic coefficient k and negligible mass are attached one to P and 
the other to the origin O. (See Figure®) It is easy to verify that the active force 
acting on P is given by 

k{y + l) ( x 
\x 2 



f(x,y) 



X 2 + 1 

Clearly /|m is a vector field tangent to M, and the motion of P is described by the 
equation £„■ = /(£) on M . Here again m = s = 1 and (J = lxl. Let 

kx(y + l) 1 



T(x,y) 



y- 



x 2 + l '2" 

which is clearly defined on U, and let fl = [0,oo) x C^(M). Since deg(J r , U) = 1, 
Theorem \4-l\ shows that if any T -periodic force \h(t,£,£), A > 0, is added to f, 
then there exists an unbounded connected set T of nontrivial solution pairs for (ll.ip 
in Q whose closure in Q meets / _1 (0)nfi = {(— \/2, -1), (0, -2), (y/2, -1)}. 
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5. Applications to a class of Differential-Algebraic Equations 

Let us now consider applications to semi-explicit DAEs of the following forms 
on an open set U C R m x R s : 

( k n / x = f(x,y,x,y) + Xi)(t ) x,y ) x,y), 

[ ' \Q(x,y) = 0, 

and 

x = \t)(t,x,y,x,y), 



(5 ' 2) \a(x, y ) = o, 

where we assume that f : TU -> R m x R s and t): K x TE/ -> K m x R s are continuous 
maps, t) is T-periodic in the first variable, and g: U — > M s is C°° and such that 
&2Q{x,y) is invertible for all (x,y) £ U. 

Let us discuss briefly the notion of solution of DAEs as (|5.1[) and (J5.2D . In 
general, let £ : R x TE/ — > R m be a continuous map and let g be as in (|5.ip and 
()5.2p . Consider the equation 

x = £(t,x,y,x,y), 



(5 ' 3) ' g(x,y) = 0. 

By a solution of (I5.3[) we mean a pair of C 1 functions x: J — >• R fc and y: J — > 
R s , J an interval, with the property that x(t) = £(t,x(t),y(i),x(t),y(t)) and 
g[x(t),y(t)j = for all t € J. Differential-Algebraic Equations as (15.31) are closely 
related to second order differential equations on manifolds as described in Section 
[3] It is now our aim to clarify such a connection. 

Let (x, y) be a solution of (|5.3[) . Differentiating twice the relation g(x(£), y(t)j = 
with respect to t, we get 

d 1 Q(x(t),y(t))x(t) + d 2 g(x(t),y(t))m + &"(x(t),y(t))(x(t),y(t)) =0. 

Thus, recalling that x(t) = £ (t, x(t) , y(t) , x(t) , y(t)) , we get 

y(t) = - [9 2 g(x(i), y(t))] _1 (d lQ (x(t),y(t))£(t, x{t), y(t),x(t),y(t)) 

+ g"(x(t),y(t))(x(t),y(t)) 
For (eR,p,iie M m and ? ,uel s (hence (p, g, w, u) e TC/), let 

/ £(t,p,q,u,v) \ 

£[t ' P ' 9 ' "' V) = ^- [ftflfo «)] "* (9ig(p, q)£(t,p, q, U, v) + g"(p, g)(«, »)) J 

The following proposition shows that DAEs as (|5.3[) can be regarded as second 
order ODEs on differential manifolds as discussed in Section [3] 

Proposition 5.1. Equation (|5.3[) is equivalent to the following second order ODE 
on M = g- 1 (0): 

(5.4) £* = ?&&£,£), 

in the sense that (x(t),y(t)) is a solution of (|5.3[) if and only if £(£) = (x(i),y(t)) 
is a solution of (|5.4[) . .Here Pfp «) denotes the orthogonal projection of W n x M s 
onio Ttp tq \M. 

Proof. Clearly any solution of (15.31) that meets M satisfies also (I5.4[) . 

Let us prove the converse. According to Equation (I3.7P one has that for any 
M-valued C 2 curve t *-* £(£) one has that, for any t, 

?&)&)= r(£(t),£(i)) 
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where Vh t -. denotes the projection of the ambient space M. k = M. m x M. s onto the 
orthogonal complement (T^^M) 1 - of T^M, and r: TM — > R fc is the map given 
by 

r(Z,v) = (&'(0\t,m^) Y(Ofa,»?) € (T^M) 1 - 
for all £ € M and 77 G TjM. Let us now see what relation exists between r and £. 
Observe that, given any M -valued curve £ for which £(£) = £(t, £(£), £(£)) , applying 
Vh t \ on both sides of this equality we get 

r{m,i{t)) = ^ t )l = *&/(*. eo. £(*)) . 

whenever defined. The arbitrariness of £ shows that 

for all £ G M , r) € T ? M and t G K. 

Let now £: J C 1 -> [/ be a solution of (|5.4[) . From equation (|3.5p we get that, 
whenever defined, 

m =-p m £(t,mm) +r(m,m) 

Hence, 

ew =?«*)%£(*),&*)) +?w(«(*U(*)) = s(t, m,i(t)), 

and the assertion is proved. □ 

We observe that the equivalence of (J5.3I) and (|5.4p could be proven also using 
the explicit expression for r found in Section [3] (see formula (|3.11jl ). 

The equivalence shown by Proposition 15.11 allows us to apply the results of the 
previous sections to Equations (|5.ip and (|5.2[) by the means of the corresponding 
ODEs on M. 

Let f, g, and f) be as above. Define the maps *: T(U) -> W n x R s and T: R x 
T(j7) -> K m x M s by 

&(p,q,u,v) = 

(f(Pi 9, «,*>),- [§2fl(p, <?)] ^igO, q)f (p, q, u, v) + q"(p, q) (u, v) ^ 

and 

T(t,p,q,u,v) = 

[t)(t,p, q, u, v), - [d 2 g(s>, q)] 1 [diQ(p, q)t)(t,p, q, u, v) + g"(p, q) (u, v) 

We have that (|5.1[) and (|5.2p are equivalent, respectively, to 

(5.5) ^ = ^*(^) + A^T(£,£,0, 
and 

(5.6) ^ = A^T(i,£,0, 

where £ = (x,y). 

We now use Equations (J5.5I) and (J5.6P to derive, by the means of Theorems 14.11 
and 14.21 results concerning the sets of solution pairs for Equations (|5.ip and (|5.2p . 

Let us consider first Equation (|5.1JI . We define fo(p, (?) := f(p, q, 0, 0) and 

*o(p,?) :=*(?,«, 0,0)= (fo(p,g),-[5 2 0(p,g)] _1 9i0(p,g)fo(p,g)) • 

Observe that ^o is tangent to M in the sense that \I>o(p, q) G T(p,q)M f° r an y 
(t,p,q) G M x M . Hence, 

Pri"P( P , 9 )*o(p,g) = fo(p,?)- 
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Then from Theorem 14.11 and Remark 14.31 we get the following. 

Corollary 5.2. Let f, f), and U be as above. For any (x,y) € U, define J- : U — > 
R m x W by 

F(x,y) = (fo(x,y),g(x,y)). 

Let ft C [0, oo ) x C\ (R m x R s ) be open. Assume that deg(J r , U fl fi) is well defined 
and nonzero. Then there exists a connected set T of nontrivial solution pairs for 
(15. ip in Q whose closure in Q meets .F _1 (0) n £1 and is not compact. In particular, 
ifQ~ l (0) «s a closed subset ofR m x R s and Q = [0,oo) x C^(K m x R s ), then T is 
unbounded. 

A similar argument allows us to deduce the following consequence of Theorem 
14.21 We focus on Equation (|5.2[) . We define f)o(t,p,q) := t)(t,p, q, 0, 0) and the 
tangent vector field Wf, on M as in (|3.13[) with £ in place of (p, q). 

Since, clearly, w-p, ,x£ T^^M for any (p, a) e U, one has that 

Pr l« ; 7'(p,,) T (P'3) = w t)o(P><l)- 
Hence from Theorem 14.21 and Remark 14.31 we get 

Corollary 5.3. Let f) and q be as above. For any (x,y) € U, define J- ': U — > 
R m x R s &u 

J"(a;,j/) = (w[, (a;,y),g(a;,y)). 

Lei ft C [0, oo) x Cy(R m x R s ) 6e open. ^4ssMme i/iai deg(J r , C7 n fi) is we/Z defined 
and nonzero. Then there exists a connected set T of nontrivial solution pairs for 
(I5.2[) m Q whose closure in fl meets .F _1 (0) n fl and is not compact. In particular, 
«/fl _1 (0) «s a closed su&sei o/R m x R s and ft = [0,oo) x C^(W m x R s ), then T is 
unbounded. 

We conclude this section with a few illustrative examples. 

Example 5.4. Consider the setting of Example \4-5\ The horizontal acceleration 
experienced by P is — fc ^j-rp ■ Hence, the dynamics of Example \4-5\ is described 
also by the following DAE: 

x — ^ J.2 + ! , 

±a: 2 - y - 2 = 0. 
Let ms perturb this equation as follows 

_ ;„ (y+l)s 



p + Af)(*,x,j/,i,j/), 



I x 2 - y - 2 = 0, 

where rj : R x R 2 x R 2 — >• R is a continuous T -periodic function. Let fl = [0, oo) x 

C^(R m x R s ) and put 

as in Example \4-.5\ Since deg(J r , C/flfi) = 1^0, Corollarv \5.2\ imvlu that in fl i/iere 
exists an unbounded connected set T of nontrivial solution pairs for the perturbed 
equation whose closure meets 

•F- x (0) n n = {( - v/2, -1), (0,-2), (V2, -1)} . 

T/ie result is exactly the same as in Example \4-5\ but now we do not need the second 
component of the active force, nor we need to bother about second order differential 
equations on manifolds. Notice also that the result holds true for any choice of the 
function f). 



14 ALESSANDRO CALAMAI AND MARCO SPADINI 

Example 5.5. Let m = s = 1 and consider the following DAE (the constraint is 
that of Examvle \3.3\) : 

x = X[x + sin(i)y), 

\x 2 - y 3 - y = 0. 

Let Q = [0, oo) x C^(R m x R s ) and put 

H*,V) = I ^ / (* + sin(t)v) d* , 2^ 2 -2/ 3 -y) = (», ^x 2 -y 3 - 

Since deg(J r , t/ PI f2) = — 1, by Corollary \5.3\ we get that in O there exists an un- 
bounded connected set T of nontrivial solution pairs for this equation whose closure 
in f2 meets {(0,0)} (regarded as a a solution pair). 

Example 5.6. Let m = 1, k = 3 and U = R 3 , and consider the following DAE 
(the constraint is that of Examvle \3.5\) : 

x = x-2y + \i)(t,x,y,z,x,y,z), 
z 3 + z - x = 0, 
z - y + x 2 = 0, 

where f) : R x R 3 x R 3 — >• R is any continuous T -periodic perturbing function. Let 
Q = [0, oo) x (^(R 1 x R 2 ) and put 

J-(x, y) — (x — 2y, z 3 + z — x,z — y + x 2 ) . 

As it is readily checked J r ~ 1 (Q, 0,0) = {(0,0,0)} and deg(J r , U fl f2) = — 1. Hence, 
for any 'perturbing' function f), Corollaru \5.3\ yields an unbounded connected set 
r in Q of nontrivial solution pairs for this equation whose closure in f2 meets the 
point {(0,0,0)} (regarded as a solution pair). Notice also, as in Example \5.4\ that 
the above statement holds true for any choice of the function f) . 
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